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We derive a non-perturbative solution to the Floquet-nonequilibrium Green function (Floquet- 
NEGF) describing open quantum systems periodically driven by an external field of arbitrary 
strength of frequency. By adopting the reduced-zone scheme, we obtain expressions rendering con- 
served charge currents for any given maximum number of photons, distinguishable from other existed 
Floquet-NEGF-based expressions where, less feasible, infinite number of photons needed to be taken 
into account to ensure the conservation. To justify our derived formalism and to investigate spin- 
charge conversions by spin-orbit coupling (SOC), we consider the spin-driven setups as reciprocal 
to the electric-driven setups in S. Souma et. al., Phys. Rev. B 70, 195346 (2004) and Phys. Rev. 
Lett. 94, 106602 (2005). In our setups, pure spin currents are driven by the magnetization dy- 
namics of a precessing ferromagnetic (FM) island and then are pumped into the adjacent two- or 
four-terminal mesoscopic Aharonov-Casher (AC) ring of Rashba SOC where spin-charge conversions 
take place. Our spin-driven results show reciprocal features that excellently agree with the findings 
in the electric-driven setups mentioned above. We propose two types of symmetry operations, under 
which the AC ring Hamiltonian is invariant, to argue the relations of the pumped/converted currents 
in the leads within the same or between different pumping configurations. The symmetry arguments 
are independent of the ring width and the number of open channels in the leads, terminals, and 
precessing FM islands, In particular, net pure in-plane spin currents and pure spin currents can be 
generated in the leads for certain setups of two terminals and two precessing FM islands with the 
current magnitude and polarization direction tunable by the pumping configuration, gate voltage 
covering the two-terminal AC ring in between the FM islands. 

PACS numbers: 72.25.Dc, 03.65. Vf, 85.75.-d, 72.10.Bg 



I. INTRODUCTION 

In this section, we first give the introduction to the 
phenomenon and effects that motivate our investigation 
in Sec. II A[ An overview of the attempts of and findings 
in our study is given in Sec. II Bl where the organization 
of this paper is also provided. 



A. Spin pumping, inverse spin-Hall effect, and 
Aharonov-Casher effect without dc bias voltage 

The spin-Hall effect (SHE) is a phenomenon where lon- 
gitudinal injection of a conventional unpolarized charge 
current into a system with either extrinsic (due to impu- 
rities^) or intrinsic (due to band structure^) spin-orbit 
coupling (SOC) generates a transverse pure spin current 
in the four-terminal geometry or the corresponding spin 
accumulation along the lateral edges in the two-terminal 
geometry. While the magnitude of the pure spin cur- 
rent generated by SHE in metals and semiconductors is 
rather small and difficult to control^ the inverse spin- 
Hall effect^ (ISHE) has recently emerged as the prin- 
cipal experimental tool to detect induction of pure spin 
currents by different sources. 

In the ISHE (which can be viewed^ as the Onsager 
reciprocal phenomenon of the direct SHE), a longitudi- 



nal spin current generates a transverse charge current or 
voltage in an open circuit. Experimental examples em- 
ploying ISHE to detect pure spin current include: (i) & 
pure spin current pumped by precessing magnetization 
of a single ferromagnetic (FM) layer under ferromagnetic 
resonance (FMR) conditions with detection by injecting 
the pumped current into an adjacent normal- metal (NM), 
such as Pt, Pd, Au, and Mo, or semiconductor layerji& 
(ii) spin currents generated in nonlocal spin valves; 9 (Hi) 
a transient ballistic pure spin current injected 10 by a pair 
of laser pulses in GaAs multiple quantum wells being con- 
verted into a charge current generated by ISHE before 
the first electron- hole scattering event, thereby providing 
unambiguous evidence for the intrinsic direct and inverse 
SHE. 

The spin pumpingii by precessing magnetization is a 
phenomenon where the moving magnetization of a sin- 
gle FM layer, driven by microwave radiation under the 
FMR, emits spin current into adjacent NM layers. The 
emitted spin current is pure in the sense that it is not ac- 
companied by any net charge flux. This effect is termed 
pumping because it occurs in the absence of any dc bias 
voltage. Particularly, the detection of pure spin currents 
pumped by magnetization dynamics has become a widely 
employed technique^ to characterize the effectiveness of 
the charge-spin conversion by the SHE via measuring 
the material-specific spin-Hall angle (i.e., the ratio of 
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spin-Hall and charge conductivities). The same ISHE- 
based technique is almost exclusively used in the very re- 
cent observations of thermal spin pumping and magnon- 
phonon-mediated spin-Seebeck effect^ Also, spin pump- 
ing makes it possible to inject^ spins into semiconduc- 
tors with electrically tunable efficiency across an Ohmic 
contact, evading the notorious problem- 1 ^ of impedance 
mismatch between the FM conductor and high-resistivity 
material. 

On the theoretical side, the mechanisms for convert- 
ing pumped pure spin current into charge current, due 
to a region with intrinsic or extrinsic SOC into which the 
pumped spin current is injected, have been anal yze d in 
a number of recent studies. For example, Ref. [l5| has 
shown that both transverse and longitudinal charge cur- 
rents arc generated in the four-terminal Rashba-spin-split 
two-dimensional electron gases (2DEGs) of square shape 
which is adjacent to the FM island with precessing mag- 
netization that pumps longitudinal pure spin current into 
the 2DEG. In this scheme, the output charge current can 
be increased by increasing the strength of the Rashba 
SOC in the 2DEG. 

Furthermore, the recent alternative description^ of 
spin pumping in FM|NM multilayers, which encompasses 
both the earlier considered^ nonlocal diffusion of the 
spin accumulation at the FM|NM interface generated by 
magnetization precession and the effective field described 
by the "standard model"™ of spin pumping viewed as 
an example of adiabatic quantum pumping that is cap- 
tured by the Brouwer scattering formula^ has shown 
that spin-charge conversion does not always occur and 
that the conversion depends sensitively on the type of 
spin-orbit interactions. That is, unlike in FM|NM sys- 
tems where spin-charge conversion is driven by the ex- 
trinsic SOC and assumed to follow simple phenomeno- 
logical prediction j c oc S x j s (j c is charge current den- 
sity, S is the spin polarization direction, and j s is the 
injected spin current density), the pumped charge cur- 
rents in Rashba systems were found to deviate from this 
naive formula. 

Thus, the whole phenomenon of spin-charge conversion 
after pure spin current is injected into a system with 
SOC needs to be discussed together with the origin of 
spin currents and the type of SOC employed for the con- 
version, 1 '' Here we analyze spin current generation by 
one or two precessing FM islands and the corresponding 
spin-charge conversion in two- and four-terminal meso- 
scopic rings, adjacent to those islands and patterned in 
the 2DEG with the Rashba SOC. Unlike the spin-charge 
conversion in experimental and theoretical studies dis- 
cussed above, where electronic transport in semiclassi- 
cal nature^ the device depicted in Fig. [T] involves spin- 
sensitive quantum-interference effects caused by the dif- 
ference in the Aharonov-Casher (AC) phased— gained 
by a spin traveling around the phase- coherent ring. The 
AC effect ^Sr— in which magnetic dipoles travel around 
a tube of electric charge, can be regarded as a special 
case of a geometric phase. For typical ring sizes and 




FIG. 1: (Color online) Schematics of the (a) four- and (b) 
two-terminal device setups. A precessing ferromagnetic (FM) 
island driven by microwaves pumps pure spin current into the 
Aharonov-Casher (AC) ring (where the spin-charge conver- 
sion takes place) patterned within a two-dimensional electron 
gas in the x-y plane with the Rashba spin-orbit coupling. The 
semi-infinite normal-metal (NM) leads, in the absence of any 
dc bias voltage, are attached to the AC ring to probe the out- 
put time-averaged (over a precession period) pumped charge 

s s 

(spin-q with q G {x,y,z}) currents, I L (I L q ), Ib {I B q )y In 

s s 
{I R q ), and It {It 9 )i in the left, bottom, right, and top of the 

device, respectively. The FM island is of length Lfm, and ex- 
change splitting A. The magnetization of the FM precesses 
around the z-axis with frequency oj and cone angle Q. The 
AC ring whose lattice sites along the tangential and normal 
directions are denoted by n and m, respectively. 

strengths of Rashba SOC in InAlAs/InGaAs hcterostruc- 
tures, the AC phase acquired by a (spin) magnetic mo- 
ment moving in the presence of an electrical field is of 
Aharonov-Anandan 22 (rather than Berry) type due to 
the fact that the electron spin cannot 2 ^ adiabatically 
maintain a fixed orientation with respect to the radial 
effective (momentum-dependent and, therefore, inhomo- 
geneous) magnetic field associated with the Rashba SOC. 
In fact, the AC phase for spins traveling around the meso- 
scopic ring consists of not only the geometric phase, but 
also a dynamical phase arising from the additional spin 
precession driven by the local effective magnetic field i 2 ^ 
Accordingly, giving electrons such geometric phase 
makes it possible to manipulate the magnitude of charge 
and spin currents in AC rings due to the fact that, unlike 
usual case of intrinsically fixed phases, the ring experi- 
ments allow one to steer geometric phases in a controlled 
way through the system geometry and other various tun- 
able parameters^ For example, the destructive quantum 
interferences, controlled by the accumulated AC phase 
via tuning of the strength of the Rashba SOC (which 
depends on the applied top-gate voltag o 25 ' 26 ), cause un- 
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polarized charge current injected into the two-terminal 
AC rings to diminish 2 -^ - — [to zero^i if the ring is 
strictly one-dimensional (ID)]. Similarly, in four-terminal 
AC rings one encounters quantum-interference-controllcd 
SHE, predicted in Ref. 1 33 and extended to different types 
of SOC and ring geometry in Refs. Inland 133. where spin- 
Hall conductance can be tuned from zero to a finite value 
of the order of spin conductance quantum e/4-7r. 



B. Methodology and key results 

The goal of this study is threefold: (i) to provide 
a unified microscopic quantum transport theory based 
on the non-perturbative solution of the time-dependent 
nonequilibrium Green function (NEGF) in the Floquet 
representation— which conserves charge current at each 
level of approximation (i.e., number of microwave pho- 
tons taken into account depending on the strength of the 
driving field) for both the spin current generation by the 
magnetization dynamics and spin-charge conversion in 
the adjacent region with SOC; (ii) to understand how 
output spin and charge currents from multiterminal AC 
ring device (such as in Fig. [IJ can be controlled by the 
top-gate covering the ring, by the cone angle of precessing 
magnetization set by the input microwave power driving 
the precession, and by the setup geometry; (Hi) to ex- 
amine if the device setup in Fig. [T]can be used as a new 
playground for experiment o 21 ' 24 ' 27 ' 28 measuring charge 
currents to detect quantum interference effects involving 
AC phase in a single mesoscopic ring where multichan- 
nel effects in a typical ring of finite width act as effec- 
tive dephasing (by entangling spin and orbital degrees of 
freedom^ or averaging over orbital channels with differ- 
ent interference patterns^!), thereby randomizing inter- 
ference patterns as in conventional measurements using 
dc bias voltage.— 

The paper is organized as follows. In Sec. [H] we spec- 
ify our pumping device and the adopted Hamiltonian. 
Section Hill formulates the solution to the Floquet- NEGF 
equations. Our numerical results are discussed in Sec. lIVI 
according to the chosen parameters and units given in 
Sec. IIV Al In Sec. IIVB1 we examine both the pumped 
charge and spin currents responsible for the AC phase 
and ISHE effects and driven by spin-pumping in the ab- 
sence of any dc bias voltage, i.e., the spin-driven setups as 
the counterparts to the conventional voltage-bias driven 
(electric-driven) setups with two-terminal2^— and four- 
terminate mesoscopic AC rings of the Rashba SOC. Sec- 
tion |IVD] illustrates different pumping symmetries of the 
AC ring. We conclude in Sec. [V] 

Our key results are as follows: (i) To arrive at Eqs. 
([52]). (PI), (gl]) , and (ESI) , we solve the Floquet-NEGF 
equations and use the so-called reduced-zone scheme^ 
which guarantees conservation of charge currents for any 
given maximum number of photons, unlike other recent 
approaches based on continued-fraction solutions^ 7 - - — 
where charge conservation is ensured only in the limit 



of infinite number of photons, (ii) With Fig. [2] through 
Fig. we analyze the pumped currents in the spin-driven 
setup Fig. [T] The results are in good correspondences 
to the reciprocal electric-driven results shown in Refs. [3l| 
and l32l justifying the derived formalism herein. Detail 
examinations, based on the AC effect and ISHE, of the 
modulations of both the pumped charge and spin cur- 
rents are given. (Hi) In Sec. IIV Dl we tailor the pumping 
symmetry under which the Hamiltonian of the AC ring of 
Rashba SOC remains invariant. By performing the sym- 
metry operations on one specific pumping configuration, 
we can obtain the relations between pumped currents in 
the same or different pumping configurations (or setup 
geometry). Although we illustrate the symmetry opera- 
tions by considering only the setups of two-terminal two- 
precessing FM islands, the symmetry arguments are ap- 
plicable to the case of arbitrary number of terminals and 
FM islands as well, giving multifarious manipulations of 
the pumped currents via setup geometry. In particular, 
Fig. [H Fig. [TTJ and Fig. Ql through Fig. [T7] show that 
the pumped spin currents are pure and are of magnitude 
and polarization direction tunable by the top gate volt- 
age controlling the strength of the Rashba SOC and by 
the pumping configurations. 



II. DEVICE SETUP AND HAMILTONIAN 

Consider the spin-driven four-terminal (or four-lead) 
setup in Fig. [Ua). A ferromagnet, FM, with precession 
axis along the z direction contacts the AC ring of Rashba 
SOC in the x-y plane from the left. The FM plays the role 
of a spin-z source, pumping pure spin-z currents into the 
ring via the FM|AC-ring interface. The spin-charge con- 
version takes place in the AC ring. The pumped or con- 
verted charge current I p and spin current I p 9 are probed 
by the NM leads where currents are conserved with p = 
L, R, B, and T indicating the currents flowing through 
the left, right, bottom, and top leads and q e {x,y,z} 
standing for the pumped spin- a:, y, and z currents, re- 
spectively. 

All computed pumped spin and charge currents are 
time- averaged (over one precession period); they carry 
positive signs if the flow direction is in +x or +y direc- 
tion or minus if flow direction is in — x or —y. In the 
two-terminal setup, Fig. [ljb), we have two NMs labeled 
by p = L, R. Note that, except the number of leads, 
Fig. QJb) does not differ from Fig. [Ua), but just fur- 
ther shows the lattice structure of the device. The AC 
ring is modeled by m — 1 • • • M concentric circles of the 
same number of lattice sites, and in a circle m the lattice 
sites are indexed by n = 1 • • • N. For instance, we have 
(M,N) = (3,8) in Fig.[]Tb). The NMs and the FM are 
modeled by square lattices, while each NM is of semi- 
infinite length, and the FM is of finite length, namely, an 
island. 

The Hamiltonian of the whole device can be divided 
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into six terms, 

H (t) — Hacr + H^m + Hfm (t) 

+-HnM-ACR + #FM-ACR + ^FM-NM , (1) 

where Hacr, ^nm, and -Hfm account for the Hamilto- 
nian of the AC ring, NM, and FM, respectively. The 
term -Hnm-ACR describes the hybridization between NMs 
and AC ring, and i?FM-ACR (^fm-nm)i the hybridiza- 
tion between FM and AC ring (FM and NMs). Note that 
the time-dependent Hamiltonian originates only from the 
precessing FM, f/pM^)- Below, we express these six 
terms explicitly. 

Focus on Hacr first. As given in Ref. H2, the ring 
Hamiltonian can be written as, 



Hacr, 



E 



, n,m:CT u '«,m:o' 



<T,er'=-|\i 

(N M 
\ \ ^ n,7i+l,?n;er,o"' * + 
/ j / j Iff, a n,m;<T a- n+l,m;<T' 
n—1 m— 1 

N M-l \ 

E\ ^ m,m+l,n;a,cr' -f A 1 

/ y fr u n,m;a u, n,m+l;a' 

n—1 m— 1 / 

+H.c, 



(2) 



with n and m denoting the lattice sites along the tan- 
gential (<fr) and normal (r) directions as illustrated in 
Fig.QJb), respectively. The creation (annihilation) oper- 
ator at site (n, m) of spin a is ajj m . a (a n , m;o -). The on-site 
potential e n ' m at site (n, m) takes into account the dis- 
order and can be tuned by applying a top-gate voltage. 
In what follows, unless further specified, we will assume 
that the AC ring, NM, and FM, are all clean conductors, 
i.e., of zero on-site potentials. The hopping along the <f> 
direction, 



71,71+1, m 



1 



7o h 



(r m / a yAcb 2 
7so 

— >l 

(r m /a) A(j) 

X (a x COS 4> n ,n+l 



a y sin 4> n ,n+i), (3) 



and along the r direction, 



7 ; 



m,m+l,ra . 



7 / s + ?7so [<J y cos (f> n - a x sin <j> n ) , (4) 



consists of two terms proportional to 70 that originates 
from the kinetic energy and to 7so that results from 
the Rashba SOC, with <fr n = 2ir (n — 1) /N, (f> n ,n+i = 
(4>n + K+i) /2, A</> = <p 2 - <Ai, r m = n + (m - l)a, 
70 = h(2ma )~ , 7so = a (2a) , a being the lattice 
spacing, a q = S q 2/h being the Pauli matrices, and Hx2tt, 
the Planck constant. Being worth addressing, the Hamil- 
tonian ([2]) yields the same spin precession as obtained 
by the SU(2) non-Abelian spin-orbit gauged that ab- 
sorbs the Rashba SOC term for the U-shaped 42 ID con- 
ductor; furthermore, the above form of the concentric 



tight-binding Hamiltonian was also used to theoretically 
model the Rashba SOC in HgTc/HgCdTe quantum wells 
in Ref. H3, showing experimental observations of the AC 
effect in good agreements with the theoretical predic- 
tions, and thus strengthening the validity of the ring 
Hamiltonian Eq. @. 

The currents are probed by the un-biased NM leads 
whose Hamiltonian reads, 



P <r=t,4- (l*,lt') 



(5) 



where bffij is the creation operator and b$ a is the anni- 
hilation operator in lead p at site /i of spin a. The pure 
spin currents are pumped by the precessing FM described 
by, 



Hfm (t) 



-EE 7o4;<t^'; CT 

E E [ue i( " t+$) + yt e -*M+*) 



CT,cr'=t-l v 

A 

+ — cos Oct, 



■ E E 7(4;^ 
<r=f,4. {vy) 



(6) 



with M (t) = [sin O cos (wt + $) , sin 9 sin (ut + $) , cos O] 
giving V = sin©^ -ia y )A/A. The cj, ;(7 (^ ; <t) is the 
creation (annihilation) operator at site v in the FM of 
spin ct. The hybridizations between adjacent materials, 



H 



VVl-ACR = -70 E E E 'V<:""<<- + H.C., 
P cr=t4 (jU,n) 



#FM-ACR = -70 E E + H ' C -' 



and 



FM-NM 



7>E E E & l,*h;ii° + ll r - 

p CT=f4 



are set to be of the same strength, namely, 70. 



III. FLOQUET-NONEQUILIBRIUM GREEN 
FUNCTION APPROACH FOR PERIODICALLY 
DRIVEN OPEN QUANTUM SYSTEMS 

In the devices where spin flip or spin precession is ab- 
sent, the problem of spin pumping by magnetization dy- 
namics can be greatly simplified by mapping it onto a 
time-independent one in the frame rotating with the pre- 
cessing magnetization42r— However, the device in Fig.Q] 
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contains Rashba SOC which causes spin-up to evolve 
into to spin-down by spin precession, so that the device 
Hamiltonian transformed in the rotating frame contains 
time-dependent SOC terms. 

In the adiabatic regime u> —> 0, which is satisfied for 
pumping by magnetization dynamics since the energy of 
microwave photons hu is much smaller than other rele- 
vant energy scales j* 7 - one can employ the Brouwer scat- 
tering formula. 18 However, this is numerically very in- 
efficient since all pumped spin and charge currents in 
devices, where the precessing FM island is coupled to a 
region with SOC, are time-dependent Thus, one has 
to compute scattering matrix of the device repeatedly at 
each time step of a discrete grid covering one period of 
harmonic external potential in order to find full ac cur- 
rent vs. time dependence and then extract experimen- 
tally measured dc component. 

The relevant dc component of pumped current can be 
obtained from approaches which generalize the existing 
steady-state transport theories, such as the scattering 
matrix^ 8 .^ NEGF formalism ) 35 ' 37 ' 38 i 50 ~-^and quantum 
master equations 5 ^ with the help of the Floquet theo- 
rem's ygjjd f or periodically driven systems. While the 
equations of the Floquet-NEGF formalism we adopt here 
have been used before to study a variety of charge pump- 
ing problems in non-interactin g 37 ' 50 " — and interacting 
electron system o 35 ' 53 or the photon-assisted dc trans- 
port j 3 - 8 - the key issue is to find a solution to these equa- 
tions that can capture pumping processes at arbitrary 
strength (or frequency) of the external time-periodic po- 
tential while conserving^ charge currents at each step 
of analytic or numerical algorithm. For example, the 
often used continued fraction solution 3 - 7 - - — to Floquet- 
NEGF equations does not^ 7 - conserve charge current in 
the leads, and the key trick we employ below to ensure 
current conservation is the reduced-zone schemed 

We begin the derivation for the charge-current- 
conserved Floquet-NEGF solution by noting that the two 
fundamental objects 5 -^ of the NEGF formalism are the re- 
tarded 

G^j (t, t') = - i« (t - {[d T (t) , dt, (*')}) (7) 



and the lesser 



G< (t,f)= % - (4 (t')dx(t) 



(8) 



Green functions which describe the density of available 
quantum states and how electrons occupy those states, 
respectively. Here u is the unit step function; indices 
{X, J7} £ {n, m, /x, is, crjand creation d) or annihilation 
operators d £ |a, b, c| are used. For notational conve- 
nience, the matrix representation with indices {X, J7} will 
not be written out explicitly below. 



The essence of the Floquet-NEGF approach is to treat 
the time variable t in Eq. (JTJ) as an additional real-space 
degrees of freedom denoted by t with considering the aux- 
iliary first-quantized Hamiltonian 



h(t) = h (t) - ih 



d_ 
dt 



(9) 



Here h (i) is the first-quantization version of our actual 
or original Hamiltonian ([T]), i.e., the matrix represen- 
tation for h (t) is of elements hxj (t) given by H (t) = 
Si j d\djhzj (t). The check-hatted symbol X is used 
to remind us that X is an auxiliary variable or operator 
but not the actual one. 

The Schrodinger equation for h (t) reads, 



d T 

-at 1 



ih^-i>{i,t) =h(t)i>(t,t) 



(10) 



while keeping in mind again that only t is the real time 
variable, but t is a virtual position variable. It is straight- 
forward to prove that, by assuming the wave function 
tp (t, t) of the form tp (t, t) = A (t) B (t) in Eq.flTU]) and 
then letting t — > t, the original wave function is recovered, 



$(i=t,t)=il>(t), 



(11) 



where ip (t) obeys our original Schrodinger equation 
ihdip (t) /dt = h (t) ip (t) . Equation (fTTj) plays the fun- 
damental role in the Floquet-NEGF, since it bridges the 
two systems, the auxiliary time-independent system de- 
scribed by h (P) and our original system described by 
h(t). Accordingly, one can first solve the problems in 
the time-independent system constructed according to 
Eq. ((9]) , express physical quantities or functions in terms 
of ip (i, t) , and eventually set t — > t to obtain the corre- 
sponding physical quantities or functions for our original 
system. 

To illustrate the idea above, consider the retarded 
Green function as an example. The retarded Floquet 
Green function G r (t, t'; i, Fj corresponding to our auxil- 
iary system h (t) obeys the equation of motion (EOM), 



G r (t, t'; i, i') =S(t- t 1 ) S T (I - F) 



(12) 

where St (t — i') denotes the Dirac delta function of pe- 
riod T, Note that h (t) is time-independent; hence, 
G r \ t,t']t,t'\ depends only on the single time variable 
t — t' via the Fourier transformation 

G r (t, t'; i, I') = J°° ^_ e -^ t - t ')' h G r (E; i, F) (13) 
and can be expanded by the wave functions of the form 
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G r (E;i,t') 



E 



(*) 



n P h,m P h = — oo 



Here, the notations, identity operator 1,7]—} 
0+, {n ph ,m ph } G integers, and (■■■) nph ^ rrlph = 

I-T/2 ^n P h (*) (' ' ') VVn ph (*) are used, and the basis 

=(T)- 1/2 e-"^* (15) 

ensures the periodicity G r (E; i + IT, i' + I'T) 
with uj = 2it/T and {/,/'} G integers. The 

\[(E + it])I- ft-]" 1 ) in Eq. (HU) is evalu- 

I- J ri ph ,m ph 

ated according to the definition © of h via 



T/2 



-T/2 



^n ph (*) ^ (*) ^m ph (*) 



— ft-n p h,m P h n ph^ w ^n p h."i p h J 



(16) 



with n p h < (n p h > 0) accounting for the absorption 
(emission) processes of photons, as indicated by the sub- 



1 

(E + irj) I — h_ 



^m ph (<0 



(14) 



n p h,»Ti ph 



r 



script "ph", and 



ftn P h,f 



T/2 

-T/2 
1 /- T / 2 



T 



d ^ ph (*) & (0 ^ (*) 

d f e i(n ph -m pll )wi h ^ _ ( 17 ) 



T/2 



Note that here beside the {I, J"} degrees of freedoms, the 
extra degree of freedom, namely, photon is introduced. 
The ftn ph ,m ph exists in the {I, J} ® n p h Hilbert space, 
and thus so does G r (E] t, i') . The primary result for the 
actual (original) retarded Green function (JTJ) is obtained 
by substituting (jl"4"]l . computed by Eqs. (fl"5"j) . (jTHJ), and 
(ITTl) , into Eq. (TT3t , and then replacing t with i and t' with 
t' in the wave-function-expanded Green function, namely 
G r (*,< / ;*,* / )|f_ >ii '^t' = Gr (^ < '); nonetheless G r (M') 
can be further simplified by noting the relation that for 
l p h G integers, one has 



1 



(E + irj) I -h 



T-ph^ph+^ph 



1 



(E + Iphtu-'J + iff) ~ ft 
I 



(18) 



which can be deduced simply from the observa- 
tions that A4 = (E + ir)) I — h is a matrix of in- 
finite size, and that \(E + irf) I — h\ , and 

[(E 1 + Lhfiw + 177) — /il , evaluated according to 

(fT6"|) and p7|) are at the same position of A4, i.e., the 
same matrix element of A4] therefore, the element at 
the same matrix position of Ai^ 1 implies Eq. (TT51) . 
a manifestation of the reduced-zone scheme in which 



"ph - ' P h,™ph 



the energy E can be reduced to the zone of range 
hw. With the help of ^[(E + irj) I - h]^ 

<[(E + m pil huj + irj) —h\ 1 ) 



-m ph ,0 



'ph 



P h,"i P h 

= and 



the re-notation Z p h — > m p h in Eq. (1181) ] and change of 

^ P h - TOph, the 



variables, E' 



E 



hw and 



retarded Green function now reads, 



C'(M') = & (t,W)\ M ^ t , 

= I f°° dE c -iE(t-t')/h -m ph 

T J_ oc 2ttH ^ 

n ph ,m ph = -oo 

/OO J T?l 00 
J^f iE ' {t ~ t,)lh E e-^G^AE'), 



(E + irj) I -h 



™ P h,™ P h 



(19) 



with 



G^ mp jE) ee {[(E + l r 1 )I-h]- 1 } r 



G 



q (E + rn p hfiw); note that the prefactor S (0) in 



7 



(fT9| is due to the cancelation of m v \Jujj in the exponent, 
i.e., 6(0) = (T)' 1 J2 mph 1 = m" 1 E mpi e^ xG , 
while for physical quantities, this prefactor S (0) is ir- 
relevant; instead, it is the normalized [absence of <5 (0) in 
([TTJ]) ] Green function, 



dE 

2^h t 



-iE(t-t')/h 



e-^G^ fl {E), (20) 



n T ,- h = — oo 



that renders physical observable. One can also verify that 



J 



the expression (|20j) satisfies the EOM, 



m- -h{t) 



G T {t,t') = 5 T (t-t'), (21) 



by applying /5r/2 ^ t° both sides of the above Eq. ([2 ] 



-T/2 ' 

Comparing the EOMs (|12[) and (|2"Tj) , one clearly sees that 
the evolution of G r (i, i') is governed by the actual system 
h (t), while G r (t, t'; i, i') is by the auxiliary system h (J) . 

The lesser Green function can be obtained in the same 
manner. The auxiliary lesser Green function obeys the 
Keldysh integral equation, 



/oo /.oo pT/2 pT/2 

/ dt x dt 2 / / dhdi 2 G r (i, ti; t, h) S< (ti,t 2 ;h,h) G a (t 2 ,t'; t 2 ,i') ■ (22) 
-ooJ-oo J -T/2 J -T/2 



I -T/2 J -T/2 



r 



The time-independent h(t) allows G K (t,t';i,i') to be with G < (E;i,i') being Eq. (j2"2")l written in the energy 

expressed in terms of the single time variable t — t of the domain, 

form, 



G< (t,t' ;t,F) 



dE 
2nh t 



iE(t-t')/H G < ( E .j ^ 



(23) 



f T/2 f T/2 

G< (E; t, t') = I / dhdi^ (E; i, ij) £< (E; i u i 2 ) G a (E; t 2 ,t') 

-T/2 J —T/2 



(24) 



r 



Here G a (E;i,t') = [G r (E;i' ,i)y is the advanced 
Green function, and t< (E;i,i') = |7o| 2 ff < (E;i,i') = 
J2 P ItoI 2 ^^ (E;t,t') is the lesser self energy account- 
ing for the interactions from all probes with the bare 
(probes that are free of interacting with the environ- 
ments) lesser Green function of probe p denoted by 
g(p)< (E;{,P). The primary result for the lesser Green 

function ([8]) is obtained via G < (t, t'; i, i') L { , — 



G < (t,t r ) where G < (t,t';i,t'j is computed by the Green 
functions in the time-independent system h (t) according 
to Eqs. d23J) and Similarly, G< (t,f) can also be 

further simplified by taking advantages of relation (jT5J) 
and change of variables. For this simplification, we utilize 
the Keldysh equation in the energy domain Eq. (|24[) and 
the wave-function expansion to obtain the expression of 
G<(t,t';U'), 



6< < ^ • ^ • ^ • f o = T/„lr' E(, -'' ), * E E 

n p h,m p h = -oo fc ph ,i ph = -oo 



T 



[E + in)I -h 



« P h,fc P h v 



(E-irj)I -hi . 



T 



e -i(E'-m ph tuo)(t-t')/h 

2nh 



8 



E 



n D h,"iDh = — oo 



E 



1 



(E' — rriphHu} + irf) I — h 
1 



(£' -irf) I- ht 



(25) 



r 



where | [(J3 -irf)I - ht] 

I. J (ph,™ P h 

< [(£7 + TOphfiw — irj) I — hf\ > and change 

of variables, E' = E + m^tuo and Z ph = Z p h — "i p h 
are used. Note that the original Hamiltonian 
([5]) of the NM probes in our present case is 
time-independent so that we have the expression, 

l7o| 2 3S%p h (£) = -2»/£ ( S ) l7o| 2 Im5i P P l r ,m P h - 



,■ f (p) / p\ f(p) 

'Jriph 1-°; 1 n P h,m P h 



(of Fermi energy Ep at zero temperature as the regime 
we are interested in), 



(£), with the Fermi-Dirac distribution and 



lim 

/3-yO 



I i e (E+n ph hu>~E F )/P 



(26) 



fW )mp j£)=i| 7 o| 2 <j ( J B + ^)/-/ t (p) 



(£ -irj)I- 



(27) 



n P h,m P h 



Here the definition © yields /i^ (£) = - ift<9/<9f 
with being the first-quantized version of i?NM — 
Ylrph 1 bjfip, and note again because £Tnm in Eq. (O 
or is time-independent, one has \h^] = 

^n P h,m ph |A (p) ] nph > resulting in f^™^ (£) proportional 
to 

^n P h,mphJ i- e -i diagonal, and thus the bare lesser Green 
functions of NMs are diagonal in photon (or Floquet) 



(jT5) to ffi^mp!, (E) evaluated by Eqs. (26} and (p7|, one 
deduces, 



*(p)< 

'ri ph> m ph +Z ph 



(£) 



X(P)< 



ph _ tpW'tph 



(£7- 



'phj 



(28) 



space g^ mph {E) 



x( P )< 
■'nph,m p h.ynph 



(£) as well. More- 



over, applying the same argument we used to derive Eq. 



J 



which reflects again the reducible property (energy E 
can be reduced to the zone of range fku) that yields 
the reduced-zone scheme. Using above relation (|28|) and 
change of variables, k' ph = fc p h — m p h and rx ' h = n^— 
in Eq. (|25]l. we arrive at, 



G < (tt'-it') = — c -i(E' -m ob hu)(t-t')/h e -i(n ph +"i p h)"* e "n P h"t' 

-£ J —oo 27rfi, 

" ph ,m P h 

x y \ — - — -1 y ho\ 2 g^K {e') 

J-!, [(E' + irf)I-h\ , ., ^ Xm v ' 

fcph.'ph LV " J »ph. fc P h p 



(£!' - irf) I - /it 



(29) 



The actual lesser Green function can be obtained again so that the actual lesser Green function can be written 
via setting i — > t and P — > t' in Eq. (f2U| and noting as, 
that rriphhuj in the exponent is canceled out; we thus 
have 6(0) = (T^E^l = (T)" 1 £ mph e --P^*o, 



G< (t,t') = G< (t,t';i,i')\ 
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r ^ e -<B'(t-t')/s'r e -i<h«* v 

ph ph 



.(E' + ir))I-h 
and eventually obtain the normalized lesser Green function 

Q<(t,t') = J°° ^e-^M/^e^ £ 



(£' - ir?) J - /it 



ph' 



1 



(£7 + irj)I - h 



^ph ^ph 

Ei7oi 2 ^S Ph (^) 



(£7 - irj) I - /it 



(30) 



i P h,0 



r 



Physical quantities can be extracted from the ac- 
tual Green function (|30|). For instance, the quantum- 
statistical-averaged occupation number at time t on site 
i e {n,m,fj,} can be expressed as Tr s \Gf i (t,t)0] H/i, 
while the bond current from site i to site j reads, 



Ji-^j (t) = -Tr s 



with O = I s (O = S q ) for particle (spin S q ) occupations 
or currents and I s being the 2x2 identity matrix in 
the Pauli space; the notation Tr s stands for performing 
the trace in the Pauli space (spin Hilbert space), and 
the anticommutator {A, B} is defined as AB + BA. The 
particle (charge) current 



%%<(M) 



(31) 



M') = ^£ E 



2nh 



E F +hui/2 



dETr' 



E F -hu/2 



p "ph,™ P h 

G r (E) f (E) G a (E) /GO (E) f GO (E) 



G r (E) f( p ') (E) f (E) G a (E) f GO (£7) 

-i(n ph -m ph )wl 



and spin current 



4 s ' (*) 



4tt 



E E 



Ep+hui/2 



dETr' 



E F -hu)/2 



G r (£7) f ( p ') (#) G a (E) /GO (£7) f GO (£7) 

I 



CTq G r (E) f( p ') (E) H p ') (E) G a (E) f GO (E) 

-i(n ph -m ph )uit 



probed by lead p are obtained by summing all the bond 
currents (f3"Tj) flowing through lead p. Note here the no- 
tation Tr' performs the trace over all degrees of free- 
dom, except photon's. All the functions within the 
trace are matrices; for example, the Fermi-Dirac distri- 



(32) 



(33) 



bution is a matrix, /GO (E) = fi P p ] h (E) 

On pi ,m P h) com- 
puted by (|26p . The time-averaged currents can then 

72 
T/2 ' 



be obtained via I p =sign(p) x 2-kKT 1 J_^/ 2 dtl p (t) and 



Ip " =sign(p) x AttT- 1 f_l? /2 dtlp" (t) as 



pE F + hu/2 

Ip = sign (p) x V / d£Tr G r (£7) /( p ') (£7) f (£7) G a (£7) f GO [E) 

, JE F -hui/2 L 

- G r (E) f ( p ') (£7) G a (£7) /GO (J5) f GO (£7) 



(34) 
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and 



i-E F +hui/2 

sign (p) /_ ^ d£Tr[a 9 



Ef-Hu/2 



G r (E) /( p ') (£) f ( p ') (£) G a (£) fW (£) 



f (J5) G Q (J5) / (p) (£) f (p) (£)] } 



(35) 



r 



with sign(p) = 1 for p G {R,T} and sign(p) = —1 for 
p G {L, B}. Notice here the appearance of sign(p) is 
merely for the sign convenience, positive for right- or up- 
flowing currents, while negative for left- or down-flowing 

currents. The prefactor 2irh (An) for I p (I P q ) are adopted 

s 

so that the units of I p and I p q are the same. In other 

words, if I p measures the number of charge quanta flow- 
s' 

ing through the lead p per second, then I p q measures the 
number of spins S q quanta flowing through the lead p per 
second. The trace here now is taken over all degrees of 
freedom, including photon's. 

We emphasize that, in Eqs. {33J, ([23]), flMJl, and 

is 



(|55|) . the reduced-zone scheme^ such as (fT5| or 
adopted, so that the original integral interval [-co, go] 
over energy E is reduced to [Ep — tioj/2, Ef + tiui/2], 
and with this scheme, for any given maximum n p h, 
charge currents are conserved, namely, J2 P CO = or 
X)pSign(p) I p = 0. In principle all integers n p h should ac- 
count for transport, i.e., transitions involving any num- 
ber of photons have to be taken into account, nonetheless 
when the strength of the time-dependent field is small, 
only few photons can be absorbed or emitted by elec- 
trons near Fermi level, and thus considering transitions 
between channels of few photons are sufficient enough to 
get accurate results of currents. In our following calcu- 
lations, |riph| < 2 is chosen, since we find |n p h| < 2 and 
|^ph| < 3 do not yield significantly discernible results. 



IV. RESULTS AND DISCUSSION 



By Eqs. (|34|) and (|35p . we show and examine our nu- 
merical results for two- and four-terminal spin-pumping 
setups with the parameters and units specified in 
Sec. IIV Al In Sec. IIVB1 we first concentrate on the two- 
terminal case [Fig . QJh)] to see the counterpart physics 
shown in Ref. [3l| and then, in Sec. IIV C[ we investigate 
the four-terminal case [Fig. [TJa)] to unveil the phenom- 
ena dual to what were found in Ref. 32. Our discussions 
are restricted to the case of single precessing FM in Sees. 
HVBl and ITVCl In Sec. ITVDl we aim at building up the 
relations between probed currents in the same or differ- 
ent pumping configurations from symmetry perspective; 
the two presented symmetries yield invariant AC ring 
Hamiltonian, and the arguments on the relations based 
on the symmetries are generally capable of setups of ar- 
bitrary number of precessing FM islands and terminals. 
Nevertheless, for demonstration simplicity, below in Sec. 



IIV El where our numerical results are shown to be in line 
with the predictions given by the symmetry arguments, 
we consider only the two-terminal two-precessing-FM se- 
tups. 



A. Parameters and units 

The following parameters and units are used. All en- 
ergies are in unit of the hopping energy 70, and lengths 
are in unit of the lattice constant a. For brevity, the as- 
pect ratio 1/2 between the length of FM, Lfm> an d the 
length of AC ring, N, is adopted; for example, in Fig. 
CUb), we have L FM /N = 4/8. Also, the width of FM is 
set to be the same as the width of NM. The default val- 
ues of parameters of the precession FM (FMs), the spin 
splitting strength A = 1, precession frequency (energy) 
hui = 10~ 3 , precession cone angle 9 = 10°, and initial 
precession phase (azimuthal angle) $ = 0° are chosen. 
To compare the results of our spin-driven setups with 
the findings of the electric-driven setups in Refs. HH and 
[32l . the size of the AC ring are set similar or according to 
Refs. Hi] and H2- We refer to the ring of M = 1 as the 
strict ID ring, and M > 1 as the quasi ID ring. Note 
again that the sign convention used here is, positive for 
right- or up-moving flow and negative for left- or down- 
moving flow. No bias is applied to any probes for what 
we consider here are all spin-driven setups. The number 
of open channels M opcn in the leads is adjustable by vary- 
ing Ep; referring to Fig. 2 in Ref. [U for leads of width 
consisting of three lattice sites, one has Af opon = 1 ap- 
proximately in the interval Ep G [±3.9, ±2], M opcn = 2 
in E F G [±2,±0.5], and M open = 3 in E F G [-0.5,0.5]. 



B. Single precessing FM island attached to 
two-terminal mesoscopic AC ring 

We begin with the two-terminal case of clean AC ring 
in the spin-driven setup Fig. [1Kb). Introducing the di- 
mensionless Rashba SOC strength, 



7Q7T 



(36) 



we find in Fig. [3] for (M,N) = (1,200) and Fig. 
for (M, N) — (3, 200) where only one channel is open 



(M, 



open 



1), the pumped spin-z 



current I R * 



probed by 



the right NM lead is a quasi-periodic function of Qr; 



11 




<10"' 




(M,N) = (3,200) 
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(a) T , 



M , 




-3.0 



FIG. 2: (Color online) Pumped charge and spin currents 
as a function of the dimensionless Rashba spin-orbit cou- 
pling strength Qr in the two-terminal spin-driven setup (top 
schematics) with ring size (M, N) = (1, 200) at different Fermi 
energies (a) E F = -1.8, (b) E F = -0.8, and (c) E F = -0.1. 
The ring is in contact with two semi-infinite one-dimensional 
leads in which currents are probed. The solid vertical lines 
here indicate the current modulation nodes Q R s at which 
the complete destructive interferences for spin-z take place, 
causing Ir = 1^* = 0. Since charge currents are conserved, 
II ~ Ir is satisfied. The spin-driven results here correspond 
to the electric-driven results, Fig. 3 in Ref. l3ll . 



specifically, for M opcn = 1 , by increasing Qr the I R van- 
ishes at certain Q^s, namely, the (AC-spin-interference- 
induced) modulation nodes; this behavior is akin to the 
electric-driven setup where the charge current disappears 
at these Q* R s^l 
s 

The I R modulation originates from the fact that a 
spin-z acquires some AC phase induced by the Rashba 
SOC when passing through the AC ring, and the phase 
difference between the upper-arm and lower-arm of the 
ring depends on the Rashba SOC strength; accordingly, 
gradually varying the Rashba SOC strength modulates 
the spin-z current. The condition of M opcn — 1 is sat- 
isfied when M = 1 (Fig. O, or when M > 1 (Fig. [3J) 
with the Fermi energy Ep only crossing one subband of 
the leads. In the former (strict ID), Q* R s are indepen- 
dent of the Fermi energy Ep , while in the latter (quasi 
ID), when one tunes Ep but keeps Ep in the regime 
-Mopcn = 1, Q*rS also remains unaffected (independent of 
Ep as long as Af opcn = 1 is satisfied), mimicking again 
the electric-driven setup. Moreover, since the spin-z cur- 
rent I Sz pumped by the FM is pure, if the spin-z encoun- 



FIG. 3: (Color online) Pumped charge and spin currents 
as a function of the dimensionless Rashba spin-orbit cou- 
pling strength Qr in the two-terminal spin-driven setup (top 
schematics) at different Fermi energies (a) E F = —3.0, (b) 
E F = —2.7, and (c) E F — —2.2 that yield the number of open 
channels, M opcll = 1. The ring is of size (M, N) = (3, 200) and 
in contact with two semi-infinite probes of width consisting 
of three lattice sites. The solid vertical lines here indicate the 
current modulation nodes Q* R s. The spin-driven results here 
correspond to the electric-driven results, Fig. 5 in Ref. l3ll . 



ters a complete destructive interference, i.e., not able to 
transport through the ring, then no charge currents will 
be generated in the right NM as well, providing that no 
passage of spins with different polarizations such as spin- 
x or spin-y occur through the interface between the right 
NM and AC ring as we will address below. We refer this 
types of nodes the AC- spin-interference-induced modu- 
lation nodes where I R = and Ir vanish concurrently 
at the same Q* R s, as indicated by the solid vertical lines 
in Figs. HandEl 

Note that when a pumped spin-z enters the ring, it 
starts to precess, and there are chances for this spin to 
become spin- a; or spin-?/ when leaving the ring to the right 
NM, so that nonzero charge current Ir ^ without spin- 
z current I R * = can be detected by the right NM. For 
I R = nodes involving processes as mentioned above are 
not the AC-spin-interference-induced modulation nodes 
Q* R s as being focused here, because they originate from 
precession but not interference. Noteworthily, although 
the ID ring (M = 1) in Fig. [2]and quasi-lD ring in Fig. [3J 
are both in the Af opon = 1 regime, there is an essential dif- 
ference between them. In Fig. [21 there are no evanescent 
modes, while in Fig. [3J there are two (M — M opcn = 3 — 1) 
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FIG. 4: (Color online) Pumped charge and spin currents as 
a function of the dimensionless Rashba spin-orbit coupling 
strength Qr in the spin-driven setup same as considered in 
Fig. |3j while different Fermi energies are chosen to have the 
number of open channels, A/ opcll = 2 for (d) Ef = —1.8 and 
(c) E F = -1.0 and M opcn = 3 for (f) E F = -0.35 and (g) 
Ef = —0.1. Unlike (a), (b), and (c) with M opcn = 1 in Fig. 
121 the modulation here becomes incomplete (absence of the 
AC-spin-interference-induced modulation nodes at which one 
has Ir = = 0, namely, absence of Qrs) due to the loss of 
M p Cn = 1, while the (incomplete) quasi-periodicity can still 
be found as depicted by the solid vertical lines. The spin- 
driven results here correspond to the electric-driven results, 
Fig. 6 in Ref. H 



channels that contribute to the evanescent modes; the 
Q* R nodes in Fig. [3] are thus slightly modified from Fig. 
[2j Also, the /£* can be nonzero at these Q* R s, simply 
because the FM pumps also spin-z currents directly to 
the left lead. To see the corresponding electric-driven 
results, compare Fig. [2] here with Fig. 3 in Ref. |3l| and 
Fig. E] here with Fig. 5 in Ref. [H 

For M opcn > 1, the conducting spin states become in- 
coherent or impure due to the entanglements between 
spin and orbit degrees of freedom.— The concept of en- 
semble average or spin density matrix has to be intro- 



duced to describe the transport of interferences induced 
by the AC effect. Furthermore, when more channels 
are open, the detected spin phase is obtained by taking 
into account the transport processes within and between 
each single channels, and each transport process gives 
different AC phases yielding different interference nodes 
(places where the complete destructive interference oc- 
curs); thereby, the overall complete destructive interfer- 
ence is washed out by different transport processes, form- 
ing the "incomplete" modulations (absence of Q* R nodes); 
for example, in Fig. 0] for M opcn = 2 and M open = 3 with 
N = 200, although one can still find the quasi-periodicity 
as depicted by the solid vertical lines, Ir and I R Z in gen- 
eral do not vanish concurrently. In addition, some of 
the pumped spins can be reflected in the FM|AC-ring 
interface before entering the AC ring, so that one has 

iff > I R Z . To see the electric-driven case correspond- 
ing to Fig. [4j compare Fig. [4] here with Fig. 6 in Ref. 
3lL Note that all the above two-terminal results preserve 
the conservation of charge currents with II = Ir- 



C. Single precessing FM island attached to 
four-terminal mesoscopic AC ring 

In the four-terminal setup Fig. Ufa), both ISHE 
and AC effects emerge, giving the inverse quantum- 
interference-controlled SHE. As shown in Fig. [5] with 
the ring size (M,N) — (1,100), four semi-infinite ID 
probes, and Ep = —0.05, the transverse currents obey 
Ig z = —I T Z and Ib = It for all Qr, signifying the ISHE. 
Figure [5] also shows the longitudinal currents in this four- 
terminal and the corresponding two-terminal setups. In 



the two-terminal setup, again, because of M opon = 1, 

s 

the modulation nodes Q* R s where Ir and I R Z vanish are 
found, rendering the complete interference modulation. 
In the four-terminal setup, at these Q* R s, although I R 
and I R vanish as well, while the longitudinal currents, 
II and do not vanish, and the inequality II ^ Ir 
shows up due to the presence of the top and bottom leads 
that break the longitudinal current conservation. Inter- 
estingly, the ISHE-induced Hall charge currents Ib = It 
disappear at these Q* R s, which demonstrates again the 
quasi-periodicity of the modulation and is dual (satis- 
fies the Onsager relation) to what was depicted for the 
SHE-induced Hall spin currents_[in the form of spin-Hall 
conductance) in Fig. 2 of Ref. [32|. Also note that trans- 
verse currents Ib and It decrease as increasing Q R due to 
the reflections in the interfaces between the AC ring and 
the top or bottom leads, a manifestation of the lattice 
Hamiltonian mismatch. 

To see how the width of the ring M affects the inter- 
ference, we consider the AC ring with fixed N — 100 in 
contact with four semi-infinite ID leads such that only 
one channel is available for transport in each lead, i.e., 
-Mopen = 1- Figure [5] indicates that the modulation fre- 
quency of the Hall currents Is = It for M = 2 is al- 
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2-probe 4-probe 




FIG. 5: (Color online) Pumped charge and spin currents 
as a function of the dimensionless Rashba spin-orbit cou- 
pling strength Qr in the two-terminal (top panel) and four- 
terminal (middle and bottom panels) spin-driven setups (see 
the schematics) in which the same strict one-dimensional 
ring of size (M,N) = (1,100) is considered. Each of the 
probes is one-dimensional. The complete modulation nodes 
Q^s that characterize the quasi-periodicity emerge through 
Ir = I ft* = in both two-and four-terminal setups. In the 
four-terminal setup, the transverse currents with Is = It and 
Ig" + Jy z = for all Qr reflect the existence of the inverse 
spin-Hall effect, and the quasi-periodicity can also be identi- 
fied via Ib = It — at the same Q* R s. Note that the Onsager 
relation is preserved if we compare with the finding in the 
electric-driven setup, Fig. 2 in Ref. l32l 




H 0.0 0.5 1.0 1.5 2.0 2.5 




FIG. 6: (Color online) Pumped transverse charge currents 
It = Ib at different ring width M as a function of the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the four- 
terminal spin-driven setup (top schematics) with four semi- 
infinite one-dimensional probes. Refer to Fig. 3 in Ref. [13 
for the reciprocal Onsager (electric-driven) results. 



addition, the presence of the week disorder plays merely 
the role to reduce the amplitudes of the modulation as 
also addressed in Ref. [32| for the corresponding electric- 
driven setup (compare Fig. [7J here with Fig. 4 in Ref. 

S3). 



D. Symmetry operations relating pumped spin and 
charge currents 



most double to that for M = 1, because in M = 2, 
one additional transport ring path appears. For larger 
width, the oscillations of Hall currents become vague 
since the multiple intertwined ID ring paths smear out 
the periodic behavior of the currents or average over the 
AC phase; nevertheless, the complete quasi-periodicity 
(Ib = It = at current modulation nodes Q*rs) is 
protected by the M ope n = 1 condition. The reciprocal 
features (for the corresponding electric-driven setup) are 
shown in Fig. 3 of Ref. S3). Note that the ISHE emerg- 
ing through Ig* = — and Ib — It is still preserved 
robustly against the ring width. 

Interestingly, the ISHE remains unaffected even in the 
weak disorder regime. Figure [7J plots the probed currents 
with different (weak) disorder strength W modeled by 
the random on-site potentials of the ring, namely, e n,m £ 
[-W/2, W/2}. The modulations of 1%" = and I B = 

It show that ISHE is robust against week disorder. In 



We extend our study to the case of multiple pro- 
cessing FM islands and examine the relations between 
the pumped currents. Consider the two-terminal setup 
Fig. [lib) with an additional FM island inserted between 
the ring and the right NM (as the schematics shown in 
Fig. ITU)). Let Qr (Ql) be the precession cone angle and 
®R be the initial precession phase of the right (left) 
FM. For M open = 1 at the condition Qr * Q R , we find 
that the spin-z currents probed by the left (right) lead 
remain almost constant when varying Qr (Ql) and/or 
®R bi other words, the left FM does not communi- 

cate with the right FM due to the complete destructive 
interference. Contrarily, in Fig. [51 with Ep = —1.8, 
(M,N) = (1,200) ring, two (left and right) ID leads, 
fixed Ql — 10° (indicated by the dash line) and $l = 0° 
in the left FM, at Qr = 5, i.e., the condition of the 

complete destructive interference is off, we see that the 

s s 

pumped currents probed by the left lead, I L , 1^" 1^" , 
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(M t N) = (l,XQQ) 




(a) Charge current l L 



.0 0.5 



1.0 1.5 2.0 2.5 



FIG. 7: (Color online) Pumped transverse (a) charge and 
(b) spin-z currents versus the dimensionless Rashba spin-orbit 
coupling strength Qr. The inverse spin-Hall effect (shown 
through Ib = It and Ig + 1^," = 0) survives in the weak 
disorder regime with different disorder strengths W of the ring 
in the four-terminal (each of the probes is one-dimensional) 
spin-driven setup (top schematics). The reciprocal Onsager 
(electric-driven) result is shown in Fig. 4 of Ref. 



and ifi" , significantly changes with Qr. For $^ depen- 
dence, noteworthily, we see that II and if " are not as 

sensitive to as I L " and I L V (for example, compare 
subfigures of Fig. [5] at & R ; » 135°). 

To establish a systematic analysis on the relations be- 
tween pumped currents, we inspect the device from the 
symmetry perspective. Recall that the Rashba SOC orig- 
inates from the structural inversion asymmetry^ mean- 
ing that the AC ring Hamiltonian, Eq. @ defined by Eq. 
d2]) and Eq. Q, does not remain the same by invert- 
ing the ring once. This one-time inversion asymmetry, 
however, leads to the conjecture that if one can somehow 
perform some inversion-like operations twice, then -Hacr 
might be invariant. Indeed, at least two types of symme- 




010 45 90 135 
Right Ferromagnet & R (°) 



FIG. 8: (Color online) Pumped (a) charge current and 
(b) spin- a;, (c) spin-i/, and (d) spin-z currents, in the two- 
terminal two-precessing-ferromagnet ring device as shown in 
the schematics of Fig. I10I probed by the left lead as functions 
of precession cone angle Qr and initial precession phase (az- 
imuthal angle) $_r of the right ferromagnet. The ring is of 
size (M, N) = (1, 200) and of the dimensionless Rashba spin- 
orbit coupling strength Qr = 5. The left ferromagnet is of 
Ql = 10° (represented by the dash lines) and 4>£, = 0°. The 
left and right leads are one-dimensional, i.e., the number of 
open channel M opcn = 1 with Fermi-energy Ef = —1.8. 



try operations can render invariant -Hacr- These two 
operations are illustrated in Fig. |H1 We refer the first op- 
eration as symmetry operation A (abbreviated as SOA), 
and the second as symmetry operation B (abbreviated 
as SOB). In SOA, we first invert the system with respect 
to the +x axis (x-inversion) and then invert again with 
respect to the +y axis (y-inversion); note that each in- 
version gives a A</> — > —Acj) and a a z —> —a z . After these 
two inversions, as shown in Fig. |H] we have, <p ~ ► 77 + <f>i 
A<fi -4 Acj), and 

(a x ,a y ,a z ) ->■ (-a x , -a y , a z ), for SOA (37) 

such that Eqs. ^ and (UJ) remain unaltered, conced- 
ing invariant i?ACR- For SOB, we first perform the x- 
inversion and then the replacement (a x , a y ) — > — (a x ,a y ); 
note that in SOB, the a y undergoes a y — > —a y —> cr y 
(a y — > —a y due to the inversion and then — <j y — > cr y 
due to the replacement). We thus get (refer to Fig. 0, 
cp -> 2ir - (j>, Acj) -» — A(f>, and 

(<J x ,a y ,a z ) -> {-(r x ,a y ,-a z ), for SOB (38) 

so that Eqs. © and ((H) are unchanged to yield invariant 
Hacr as well. 

For NMs, obviously, after SOA or SOB, the _ff NM re- 
mains the same, because all NMs are of the same spin- 
independent structural-inversion-invariant Hamiltonian. 
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Symmetry Operation A: SO A 











[ ~ Jt 












f? 

-cr_ 




^7 



Symmetry Operation B: SOB 




FIG. 9: (Color online) Symmetry operations A (SOA) and B 
(SOB). In SOA, the inversion with respect to the +x axis (as 
depicted by the inset, a;-inversion) is first performed to the 
system represented by the orange/gray arrow lying on the 
x-y plane and then with respect to +y axis (as depicted by 
the inset, y- inversion). In SOB, a;-inversion is first performed 
and then the replacement (a x ,a y ) — » —(<T x ,cr y ). The succes- 
sive figures show how the spin (a x , a y , <j z ) and polar angle <f> 
change after each inverting or replacing. 



Also, all the hybridizations (characterized by the same 
spin- independent hopping -70), Pnm-Acr, #fm-Acr, 
and Pfm-nm are SOA- and SOB-invariant. The only por- 
tion in the total Hamiltonian that might not be able to re- 
cover to its original form is the time-dependent Hamilto- 
nian i/pM {t), because under SOA or SOB the directions 
of the precession axis can vary. However, note that since 
what we are interested in is the time-averaged currents, it 
is the relative initial precessing phase $l — ®r that is rel- 
evant to these average currents, while the $l — $_r does 
not change under SOA nor SOB, because SOA or SOB 
are applied to the whole system (i.e., to all precessing 
FMs). Moreover, any operations or transformations will 
transfer one pumping configuration either to the same 
configuration or to another configuration; in the former, 
the symmetry argument will relate the probed pumped 
currents within a single configuration, whereas in the 
later, relate the probed currents between two different 
configurations; this will become more clear in the next 
section. Without loss of generality, in what follows, we 
choose systems originally at $l = <f>j? = 0° to illustrate 
how SOA or SOB helps construct the relations between 
different probed currents and verify these relations by 
inspecting our numerical results. 



E. Symmetry arguments applied to two precessing 
FM islands with two-terminal mesoscopic AC ring in 
between 



For the purpose of demonstration, we choose to con- 
sider here the two-terminal two-precessing-FM (left and 
right FMs adjacent to the left and to the right of the 
ring, respectively) setups, while one can apply the ar- 
gument presented below also to the ring devices consist 
of arbitrary number of terminals and precessing FM is- 
lands. In addition, since the definitions of SOA and SOB 
have nothing to do with the ring width M, ring length N, 
number of open channels M open , and Ep, the symmetry 
argument is valid for any M, N, M opcn , and Ep. Here, 
we choose (M, N) = (3, 200) and set Ep = -I.870 giving 
-^open = 2 to exemplify the symmetry operations. We 
use the notation convention A-B to describe the pump- 
ing configuration, with A accounting for the left precess- 
ing FM and B for the right precessing FM. Here, with 
{A, B} e {P , P0}, P (P @ ) stands for the FM that is 
of precession axis along +z (—z) axis and of precession 
cone angle O. For example, the schematics in Fig. [TUlis 
notated as Pio°-Plo°j m Fig. [T2]as Pio°--Pio°, an d in Fig. 
EI as P 90 o-P 90 o. 

Focus on SOA first. Consider P 10 °-Pio°, Fig. [TO] By 
applying SOA to Pio°-Pio°, the first step (a;-inversion) 
generates Pgo° -Pso° 1 while the second step (y- inversion) 
yields the swap L (left) O- R (right) and turns p8o°-P8o° 
into Pio°-Plo°j i-e., the original pumping configuration is 
recovered. As a result, in Pio°-Pio° we have, due to the 



s s 

y-inversion involved in SOA, Ijj = —Ir (or I L q = —Ir 
for all q components before any operations on spins), 
which then incorporated with the replacement (|37[) turns 



jSq 
1 L 



into (Ir 



' 1 L > i L I 



R ' X R 1 ± R )i m 



line with our numerical result, Fig. [TO] In Fig. [TT] 

(Pio°-Pio°)j by employing the same argument based on 
SOA, we obtain again the relations II = —Ir. and 
(/£•,/£", if") = (4*,4\-I%<). Being noteworthy, in 
the pumping configuration A- A such as Figs. [TOlandfTTl 
the probed charge currents vanish Ii = In = due to 
the left-right transmission symmetry^ resulting in pure 
spin currents in the NMs; this absence of charge currents 
can also be obtained by noting that the current conser- 
vation II = Iji and the symmetry argument that gives 
II = —Ir have to be satisfied simultaneously. 

On the other hand, for the left-right transmission 
asymmetric cases such as Figs. [12] (Pio°-Pio°) an d 
1131 (Pio°-Pio°)7 Ij-, an d Ir in general can be nonzero. 
Similarly, performing SOA on Pio°-Pio° (Pio°--fio o )j 
the ^-inversion renders p8o°-Pso° (-^80° -iW)) an d then 
the proceeding y-inversion gives Pio°-P.o° (Pio°--Pio°); 
hence, SOA transfers Pio°-Pio o (Pio°-Pio°) to the differ- 
ent pumping configuration Pio°-p.o° (Pio°-Pio°)- There- 
fore, the current Il-r m Pio°-Pio° (P.0 --P10 ) equals 
to —Ir.l m Pio°-p.o° (Pio°-Ao°)- Again, the rela- 
tions above for charge currents together with the re- 
placement (|37p make (I L X 



' l > ^l" > if' ) i n Pw° -P10 (-P0 
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FIG. 10: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin-x and spin-y currents probed by the left 
and right leads of finite width consisting of three lattice points 
versus the dimensionless Rashba spin-orbit coupling strength 
Qr in the two-terminal two-precessing-FM setup Pio°-Pio° 
(top schematics), i.e., the left ferromagnet and right ferro- 
magnet are of precession cone angle 10° and precession axes 
both along +z direction. 



FIG. 11: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin-x and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup Pio°-Pio° (top schematics) same as the 
one considered in Fig. [10] but with precession axes along — z 
direction for both the left and right ferromagnets. 



Pioo) identical to (if " , j|" , -if ") in Pio°-Pio° (Pio°- 
Fio°) and (if*, if", if*) in Pi °-Pio° (Pio°-Pio°) iden- 
tical to (if-, if", -if') in P 10 o-P w ° (P 10 o-P 10 o). All 
above features are again in line with our numerical re- 
sults, Fig. [12] and Fig. [13] It should be noted here that 
= 10° in our numerical calculation is chosen merely for 
the illustrations of symmetry operations. The symmetry 
argument presented above in fact is applicable for any 
cone angle 9 and even for the case of @l ^ Qr, i.e., the 
precessing cone angles for left (9^) and right (9^) FMs 
are different. Particularly, at 9 = 9l = 9,r = 90°, all 
the relations based on SGA shown above are preserved 
as well (see Figs. [14] [15j [TBI and [T7]) , while since the 
pumping configuration for 9 = 90° (precession within 
the two-dimensional x-y plane) is of higher geometrical 
symmetry than 9 = 10°, the probed currents are of ad- 



ditional relations as demonstrated below. 

Following the same procedure presented above, one can 
also apply SGB to any A-B configuration to relate probed 
currents in a single pumping configuration (if SOB does 
not generate another pumping configuration) or to relate 
probed currents between different pumping configura- 
tions (if SOB generates another pumping configuration). 
Here, we choose to focus on the special case with 9 = 90°. 
Unlike the case of 9 — 10° where we have no relations of 
the pumped currents between the two different pumping 
configurations, Pio°-Pio° and Pio°--Pio°7 at 9 = 90° the 
pumped currents in Pgo°-Pgo° can relate to the pumped 
currents in Pgo°-Pgoo. Applying SOB to Pgoo-Pgo° (Fig. 
E]) yields P 90 °-Pgo° (Fig. [p]) so that I L{R) in Pg o-Pg ° 
is equal to Il(r) i n Pso°-Pso° ■ This relation, again, in- 
corporated with the replacement ([55]) leads to the spin 

current (I jj r \, 1 umi I URV m ^o°-P90° equal to the spin 
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FIG. 12: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin- a; and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup Pio°--Pio° (top schematics) same as the 
one considered in Fig. [10] but with precession axis along +z 
(—z) direction for the left (right) ferromagnet. 



FIG. 13: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin-x and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup Pio°-Pio° (top schematics) same as the 
one considered in Fig. 1101 but with precession axis along — z 
(+z) direction for the left (right) ferromagnet. 



current (— i^fm, Zw*, — I L f R \) in Pgo°-Pgo°. The above 
predictions agree with our numerical results, Figs. I14l and 
[Pol Being worth noting, although Pi °-Pio°, Pio°-Pio°, 
Pw-P&o ) and -P90°--P90° all generate net pure in-plane 
(x-y plane) spin currents, i.e., Il+Ir = iff + Ir" = as 
predicted by SOA, the pumped spin currents at O = 90° 
are one to two order larger than those at 9 = 10° (com- 
pare Figs. [TUl and [TT1 with Figs. [T4l and [Pol) . Similar en- 
hancement of the pumped spin currents by the cone angle 
can also be found by comparing Figs. [H](Pio°-Ao o ) and 

m(Pio°--Pio°) with Fi gs- UB(P90°-ho°) and[T7](p90°- 
Pgo°)- The additional (beside what were obtained by 
SOA) relations between Pgo°-p9o° and Pgo°-P9o° can be 
obtained by performing SOB. Applying SOB to Pgo°-Pgo° 
(Fig. [161), we arrive at Pg <,-P 90 ° (Fig. [P?]) . with Il(r) 
in P900-P900 equal to Il(r) in P900-P90 , and then, by 
replacement (J55J), (if fa, I l \r)i lf(R)) in ^90°-Ao° equal 



to ( ^L(R) > I lJr) ' ~ I L(R)) in ^900-^90° , in line with Figs. 
[TBI and [TTJ We note that for all 8 = 90° pumping con- 
figurations, the pumped spin currents are pure (namely, 
li, = Ir = Q); this again can be achieved by consider- 
ing the current conservation together with the symmetry 
argument based on SOA and SOB. We emphasize that 
our results here show that the current polarization direc- 
tion can be tuned by the top gate voltage governing Qr, 
and the magnitude of the pumped currents can be con- 
trolled by the precession cone angle O, offering amenable 
manipulations on the output currents from the proposed 
devise. 



V. CONCLUSION 

In conclusion, by introducing the auxiliary system 
where the time domain is treated effectively as an ad- 
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FIG. 14: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin-x and spin-y currents probed by the left 
and right leads of finite width consisting of three lattice points 
versus the dimensionless Rashba spin-orbit coupling strength 
Qr in the two-terminal two-precessing-FM setup P900-P90 
(top schematics), i.e., both the left and right ferromagnets 
are of precession cone angle 90° and precession axis along +z 
direction. 



ditional real-space degree of freedom, we offer a plain 
access to the Floquet-NEGF formalism capable of deal- 
ing with time-periodic dynamic problems in full quantum 
approach. Particularly, by further adopting the reduced- 
zone scheme^ we derived expressions (f52j) , 
and (|35p in which charge currents are conserved for any 
given maximum number of photons. With the help of 
Eqs. and (pp|) . we reveal the physics attributed to 

the AC effect by considering the ring of Rashba SOC in 
spin-driven four- [Fig. HJa)] and two-terminal [Fig. [TJb)] 
setups. 

When the number of open channel in the leads is 
1, as a consequence of the AC effect, 



open 



one, M, 

the complete AC-spin-interference-induced modulation 
nodes characterized by I R = I R 3 = at certain Rashba 
SOC strengths Q* R s (where the complete destructive in- 



FIG. 15: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin- x and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup P900--P90 (top schematics) same as the 
one considered in Fig. 1141 but with precession axis along — z 
direction for both left and right ferromagnets. 



terferences occur) are found to be independent of the 
Fermi energy Ep in both two- and four-terminal cases 
(Figs. H and E]). 

In the four-terminal setup, the interference modulation 
is also characterized by Is = It = at the corresponding 
two-terminal Q* R s (i.e., Fig. El top panel Ir = I R = 
and bottom panel Is = It = vanish at the same Q* R s). 
Increasing the number of open channels by tuning the 
Fermi-energy to reach M opcn > 1 regime in quasi ID 
(of finite width) rings destroys the completeness of the 
modulation, i.e., absence of Q* R (Fig. [4j. Nevertheless, 
in the four-terminal case, we find that the ISHE iden- 
tified by I B = It and if* + j| a =0 (Fig. El bottom 
panel) is robust against the ring width (Fig. EJ) and weak 
disorder (Fig. [7]), and therefore, the proposed device of- 
fers a durable electrical means to measure the pure spin 
currents pumped by the precessing FM islands using the 
inverse quantum-interference-controlled SHE in the AC 
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FIG. 16: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin- a; and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup P900--P90 (top schematics) same as the 
one considered in Fig. [14] but with precession axis along +z 
(—z) direction for the left (right) ferromagnet. 



FIG. 17: (Color online) Pumped (a) charge and spin-z cur- 
rents and (b) spin- x and spin-y currents versus the dimen- 
sionless Rashba spin-orbit coupling strength Qr in the two- 
precessing-FM setup P900-P90 (top schematics) same as the 
one considered in Fig. 1141 but with precession axis along — z 
(+z) direction for the left (right) ferromagnet. 



rings. The above features based on our spin-driven se- 
tups reciprocally well correspond to the findings in the 
electric-driven setup, Refs. |3l] and HU, supporting our 
derived formalism. 

In addition to single-precessing-FM setup (Fig. [TJ, 
multiple-precessing-FM setup is studied. In the two- 
terminal two-precessing-FM setup where the ring is in 
contact with two (left and right) precessing FM islands, 
we find that the currents probed by the left (right) lead 
are independent of 0# and &r (Oj, and $i) of the right 
(left) FM under the condition of complete destructive 
interferences. In other words, the complete destructive 
interference blocks out the relation between the left por- 
tion (left FM and left lead) and the right portion (right 
FM and right lead) of our device, while this relation re- 
vives when the condition of the destructive interference 
is suppressed (Fig. [BJ). 

We also identified two symmetry operations, SOA and 



SOB (Fig. [9]), to examine the relations between currents 
in the same pumping configurations or different configu- 
rations. Performing SOA or SOB on an arbitrary pump- 
ing configuration together with the fact that charge cur- 
rents should be conserved, one can first relate the charge 
currents either in the same or in different pumping con- 
figurations, and then using Eq. (j37|) for SOA or Eq. (|38|) 
for SOB, one can further obtain the relations between 
spin currents. We choose to exemplify how the above 
procedure works by considering the two-terminal two- 
precessing-FM setup with precession cone angles 0£ = 
Q R = 10° (Figs. Q! QH [H andEP and L = <d B = 90° 
(Figs. [H [323EEH and OH. The relations predicted by 
SOA and SOB consist with our numerical results. Espe- 
cially, the net pure in-plane spin currents (for x-y plane 
with I L + I R = if' + 1^ = in Figs. HQl ED |H and [HI 
and for y-z plane with Il + Ir = if,* +Ir x = in Figs. [TBI 
and [PT)) can be achieved, and for all Ol = Qr = 90° 
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pumping configurations, the pumped spin currents are 
pure, namely, II = Jr = 0. Therefore, with employ- 
ing the spin-pumping device proposed here, the pumped 
currents can be controlled with their magnitudes and po- 
larization directions tunable via the pumping configura- 
tions (including the precessing cone angle) and the ap- 
plied top-gate voltage that varies Qr), giving potential 
applications in spintronics-based industry. 
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